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We propose a definitive test of whether plates involved in Casimir experiments should be modeled 
with ballistic or diffusive electrons-a prominent controversy highlighted by a number of conflicting 
experiments. The unambiguous test we propose is a measurement of the Casimir force between a 
disordered quasi-2D metallic plate and a three-dimensional metallic system at low temperatures, 
in which disorder-induced weak localization effects modify the well-known Drude result in an ex¬ 
perimentally tunable way. We calculate the weak localization correction to the Casimir force as a 
function of magnetic field and temperature and demonstrate that the quantum interference sup¬ 
pression of the Casimir force is a strong, observable effect. The coexistence of weak localization 
suppression in electronic transport and Casimir pressure would lend credence to the Drude the¬ 
ory of the Casimir effect, while the lack of such correlation would indicate a fundamental problem 
with the existing theory. We also study mesoscopic disorder fluctuations in the Casimir effect and 
estimate the width of the distribution of Casmir energies due to disorder fluctuations. 


I. INTRODUCTION 


The Casimir effect [l| is an experimentally accessible 
phenomenon which in most physical systems is theoreti¬ 
cally calculated by modeling metallic plates with one of 
two models: the Drude or plasma model. These simpli¬ 
fied models describe the linear response of electrons in 
the plates to an electromagnetic field. While the Drude 
model describes diffusive electrons subject to a random 
disorder potential, the plasma model describes ballistic 
electrons unhindered by disorder. These two models typ¬ 
ically provide similar predictions of the Casimir force as 
a function of plate separation, with the plasma model 
predicting a slightly stronger attraction than the Drude 
model. 

Quantative results from many experiments seem 
to favor the plasma model over a naive Drude model - 
in some ways, arguably, the more physical of the two. 
Many experiments attempt to account for the effect of 
electrostatic patch potentials in the plates, expecting the 



FIG. 1. (Color online) The geometry typically used in exper¬ 
imental measurements of the Casimir force is a gold coated 
sphere above a planar plate. Here we show the sphere sus¬ 
pended from a cantilever. We consider a lower plate of very 
thin metal with a weak applied perpendicular magnetic field. 


effect to be relevant for agreement with one model or 
the other. Several of these [1, Q find that the correc¬ 
tion due to patches would make agreement with Drude 
worse while others see agreement with the Drude 
model once the effect of patches is minimized. There is 
recent theoretical work to account for the contribution of 
patch potentials in several experiments Q , and while the 
results seem to weaken the case for the plasma model, 
the authors are careful to avoid claims that a definitive 
model for metallic plates has been found. 

In this work, we provide a new way to experimentally 
test the validity of the diffusive electron model by tuning 
an external magnetic field (or temperature from a less 
practical standpoint) in a Casimir system with a two- 
dimensional plate. The proposed experiment would be 
the typical experiment seen in Fig. [1] where the plate 
would be quasi-two-dimensional. We find a dramatic 
change in the Casimir effect between Drude model plates 
due to weak localization, shown in Fig. [21 that is just not 
seen with the plasma model. 

Weak localization (WL) is a well known and greatly 
studied effect [8l-[l^. most easily observed in low¬ 
dimensional disordered systems at low temperatures 
where quantum interference logarithmically decreases the 
conductivity of a sample with decreasing temperature. 
This is most easily understood via the simple Einstein 
relation, a = 2e^vD, where v is the electronic density of 
states per spin in the material and D is the diffusion con¬ 
stant. Weak localization provides a quantum correction 
to the diffusion constant, D —>■ D -\- SD, that is strongly 
dependent on both temperature and an applied magnetic 
field at very low temperatures. 

A fundamental assumption of the WL effect is that 
electronic motion is diffusive in nature, and its contribu- 



2 



FIG. 2. (Color online) The dependence of the Casimir pres¬ 
sure on the applied magnetic field between two disordered 
plates (one 3D and one 2D) at a separation of a = 250 nm. 
The 2D plate is described by the Drude model with the weak 
localization correction. The force is normalized by the ideal 
conductor result and is plotted for three temperatures-3, 1, 
and 0.1 K, from top to bottom. The Casimir pressure is nor- 
malized by the ideal result, Pq = — 5^^. The inset shows the 
conductivity of the 2D plate with WL correction as a function 
of the applied magnetic field, normalized by the uncorrected 
Drude conductivity, at the same three temperatures. 


tion to conductivity is calculated as a correction to the 
Drude model. Therefore, any impact found on the na¬ 
ture of the Casimir effect due to WL would apply only 
to a diffusive model of metallic plates and not a ballistic 
model; a sensitive experimental test of the effects of WL 
on the Casimir effect would provide a clear indication of 
whether a diffusive picture of electronic motion correctly 
describes the physics of the electrons in the experiment. 

The theory behind the use of a diffusive models relies 
upon performing an average over all possible realization 
of a disorder potential in the material. However, if this 
disorder average is done at the level of linear response 
instead of on the Casimir energy itself, then all effects 
from, e.g., the nonuniform nature of physical disorder re¬ 
alizations are neglected. While exact calculation of these 
neglected effects is impossible, it is possible to estimate 
whether ignoring them nonetheless gives a valid approx¬ 
imation to the Casimir energy, even if doing so causes 
discrepancies in other quantities, such as the entropy. In¬ 
deed, an important issue regarding the use of the Drude 


model in the Lifshitz formula-a typical method of calcu¬ 
lating the Casimir effect-is its apparent failure to satisfy 
the Nernst heat theorem [ij; 113 ] ) which there has been 
rigorous theoretical work to resolve [mil- Performing 
the disorder averaging at the level of the Casimir en¬ 
ergy should automatically resolve this apparent violation 
since exact treatment of any specific disorder potential 
must satisfy the Nernst theorem. 

In this paper, we first provide an introduction to the 
relevant aspects of the Casimir effect and our calcula- 
tional methods in Sec. im In Sec. Ell we then examine 
the correction to the Drude model that gives the weak 
localization effect, and in Sec. m we discuss the effect 
that this correction has on the Casimir force and provide 
estimates for the size of the effect. In Sec. E we justify 
our use of the Drude model when considering disorder in 
spite of concerns regarding the disorder averaging proce¬ 
dure by examining the effect that fluctuations in disor¬ 
der realizations would have on the Casimir energy and 
the ability to distinguish the Drude and plasma models. 
Finally, in Sec. I VII we summarize our results. 

II. CASIMIR EFFECT 

To explore the effects of WL on the Casimir effect, we 
consider a system consisting of two flat parallel plates: 
one thick plate described by the Drude model and one 
two-dimensional plate described by the Drude model 
with an additional term giving the weak localization ef¬ 
fect. With an experimental setup in the typical plate- 
sphere geometry, as shown in Fig. [TJ the gold layer on 
the sphere is thick enough to be most accurately de¬ 
scribed as a three dimensional material. Since the effect 
of weak localization in 3D materials is much weaker than 
in 2D films, we only consider the WL effect in the 2D 
plate. In addition to this system of primary interest, we 
also consider the Casimir pressure between two plasma 
plates and between two Drude plates without weak lo¬ 
calization for points of comparison. The latter of these 
will also give the expected behavior of the system in¬ 
cluding the 2D plate with WL correction at sufficiently 
high magnetic fields to completely suppress the WL ef¬ 
fect {H > 100gauss). In these cases we consider the 
same geometry, with one thick plate and a parallel 2D 
plate. For a calculation of the Casimir pressure in these 
systems we start from the well-known Lifshitz equation 
for the Casimir energy density at finite temperature [13, 
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which can be obtained by expanding the free energy of ^ denotes a sum over positive Matsubara frequencies, 
the two plate and photon system. In this expression, counting the n = 0 term with half weight. The functions 



















3 


^TM ’'"tb reflection coeflicients of plate i for 

the two polarizations of light. The subscript TM refers 
to the polarization where the magnetic field is perpendic¬ 
ular to the plane of incidence, and similarly for T E with 
the electric held. The reflection coefficients depend on 
both q and the Matsubara frequency uJn = 2TrnkBT, and 
may also depend on the other parameters of the system 
under consideration, such as the applied magnetic field 
H and additional temperature dependence. The reflec¬ 
tion coefficients can be written explicitly in terms of the 
dielectric functions of the plates, or alternatively 

in terms of the electromagnetic linear response functions 
of the plates, ni(iw„), which are related to the dielectric 
functions as 


The leading correction to the Drude result, which can¬ 
not be ignored in an accurate treatment of the two dimen¬ 
sional plate at low temperatures, comes from diagrams 
with maximally crossed impurity lines, shown on the sec¬ 
ond line of Fig. [3l which can be represented as a single 
diagram containing a cooperon. This approximation to 
the response function can be written as 

n = -h (5n, (5) 

where 511 gives the WL correction. An explicit calcula¬ 
tion of 5n in two dimensions, at low but finite tempera¬ 
ture T, and with an external ma gii etic field H, gives the 
Hikami-Larkin-Nagaoka formula^ Q, 


e(iujn) = 1 - n(iw„)/w^. 


( 2 ) 6n{iuJn-,T,H) = 


The reflection coefficients have the form 
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TM 






n(iu;n) 

^ 

n(w„) 
n(ia;„) - 2 g_L 


( 3 ) 


for two-dimensional plates, where we have defined q±^ = 
-l-a;^. For very thick three-dimensional metallic 
plates (thickness c? —>■ oo), the reflection coefficients have 
the form 




~ n(za;n) — q± + 
- n(ia;„) + q± - 


^Tb(9i i^n) 


\/q]_- n(ia;„) -|- q± 


( 4 ) 


The Casimir pressure is found from Eq. m by taking its 
derivative with respect to plate separation, a. 


^ y2 AeHDT^{T)J \4eHDTJ ’ 

( 6 ) 

T . 

where is the digamma function, D = is the dif¬ 
fusion constant, and is the electron dephasing time. 
This expression diverges logarithmically as T —> 0 with 
a sign opposite that of the Drude result, leading to a 
suppression of conductivity. It also has a very sensitive 
dependence on an applied magnetic field, becoming very 
small at moderate values of H 100 gauss) even at very 
low temperatures when the effect would be large in the 
absence of such a field. 

In two dimensions the primary dephasing mechanism 
at very low temperatures is Nyquist electron-electron 
scattering, and the dephasing time is given by [I^ 

where v = m/27rfi^ is the density of states per spin at 
the Fermi level for a two dimensional system. 



III. WEAK LOCALIZATION CORRECTION 

The basic inputs into Eq. o are the electromagnetic 
linear response functions for the two plates under con¬ 
sideration which contain all the electromagnetic prop¬ 
erties necessary for a calculation of the Casimir effect. 
The response function of a non-interacting disordered 
electron gas can be represented diagrammatically as in 
Fig. [3] The simplest approximation considers only the 
diagrams shown in the first line of the figure—those 
without impurity lines and all diagrams with impurity 
interaction ladders. In the long wavelength (i.e. local) 
limit, g —>■ 0, these terms combine to give the response 
function , from which the well 

known Drude result for DC conductivity can be found, 
^Drude _ _ . This approx¬ 

imation is valid for the one thick disordered plate we 
consider. 


IV. RESULTS 

We consider three systems of one thick plate and one 
2D plate: both plasma plates, both Drude plates, and 
most importantly, a thick Drude plate with a 2D Drude 
plate including the weak localization correction term 
given in Eq. (jG)). In all of these systems, we fix the plate 
separation at a = 250 nm and calculate the Casimir pres¬ 
sure as a function of either an externally applied mag¬ 
netic field or temperature, staying in the low tempera¬ 
ture regime where Eq. © is valid. Additionally, we set 
the elastic mean free path of the electrons in disordered 
plates to be 1 = 15 nm and the Fermi energy and effec¬ 
tive electron mass to be those of gold: ep = 5.53 eV and 
m* = l.lOmo where mg is the free electron mass M- 

We normalize all Casimir pressures by the ideal conduc- 
2 

tor result, Pq = with a = 250 nm. 
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FIG. 3. The diagrammatic expansion of 11 up to the leading 
correction to the Drude result. Solid lines represent disorder 
averaged electron Green’s functions, dashed lines represent 
interactions with the disorder potential, the shaded regions 
represent diffusons (labeled with V) or cooperons (labeled 
with C) and the circles represent current vertices. The first 
three diagrams of the third line together give the Drude result, 
while the last term gives the leading correction. The last line 
defines the renormalized vertex. 


In addition to the disagreement on the magnitude of 
the effect between Drude and plasma models, we find 
that there is qualitatively different behavior between the 
two when accounting for the effect of weak localization. 
The Casimir pressure between plasma plates has no de¬ 
pendence on the strength of the applied magnetic field, 
at least for such weak fields as we consider here, and 
only a very weak dependence on temperature in this low 
temperature regime—the change of the normalized pres¬ 
sure from 10 K to 0.1 K is a decrease of 1.7 x 10“^. In 
stark contrast, the Casimir pressure when considering a 
Drude plate with WL effects shows both a highly non¬ 
trivial dependence on even a weak applied magnetic field 
(at low temperatures), shown in Fig. [21 and also a sharp 
decrease with decreasing temperature (with no applied 
magnetic field), shown in Fig. [d] Both the temperature 
and magnetic field effects are expected when considering 
the Casimir pressure as a function of the conductivity 
of the plates. The sharp drop in the Casimir pressure 
with decreasing temperatures matches the drop in con¬ 
ductivity of the 2D plate obtained from theory, shown 
in the inset of Fig. [H and the strong dependence of the 
Casimir pressure on a weak magnetic field closely fol¬ 
lows the dependence of the conductivity of the 2D plate 
as obtained from theory, shown in the inset of Fig. |21 
and seen in magnetoresistance experiments with 2D thin 
films [nl, [III . Indeed, we find that applying a magnetic 
field of only H = 40 gauss perpendicular to the plates 
is enough to reduce the suppression of the pressure by 
approximately 40%. 

At T = 0.1 K and H = 0 gauss we find that by cor¬ 
rectly accounting for the effect of WL in the 2D plate 
the Casimir pressure is 11% less than if the 2D plate 
were described by a simple Drude model without the WL 


FIG. 4. The dependence of the Gasimir pressure on temper¬ 
ature between two Drude model plates (one 3D and one 2D) 
at a separation of a = 250 nm. The force is normalized by the 
ideal conductor result, and there is no applied magnetic field. 
The solid line is obtained from including the WL correction 
in the 2D plate, and the dashed line is the result obtained if 
the effect of WL is ignored. The inset shows the dependence 
of the conductivity of the 2D plate as a function of tempera¬ 
ture normalized by the uncorrected Drude model conductiv¬ 
ity. The solid curve is obtained from the Drude model with 
WL correction and the dashed line at 1 is for comparison to 
the uncorrected Drude model. 


correction. At this temperature and magnetic field, the 
change in the Casimir pressure from including the WL 
correction is larger in magnitude than the difference in 
the Casimir pressures predicted by the plasma model and 
naive Drude model, i.e.. 


pDrude 

_ 

pplasma 
J c 


_ pWL 
_ pDrude 


= 1.14 for T = 0.1 K, 


( 8 ) 


so the effect is large enough to be measurable for a low 
enough temperature. 

There are several ways to increase the size of the ef¬ 
fect even beyond this, the most straightforward being to 
lower the temperature even further. We also find that the 
effect can be increased by decreasing the electron mean 
free path, I, equivalent to increasing the impurity concen¬ 
tration, which can be seen by examining the dependence 
of Eq. ([5|) on the mean free path, given partially through 
the dephasing time in Eq. ©■ When considering smaller 
values of /, however, one must be sure that the impurity 
concentration is still below the limit of complete Ander¬ 
son localization, or else this model of diffusion breaks 
down. Alternatively, when considering much larger val¬ 
ues of /, which would make the effect smaller, one must 
be sure that the mean free path is much smaller than 
the sample dimension L or else the model of a disordered 
system breaks down. In an actual experimental system 
neither of these issues is likely to arise. 

Since a controlled smooth variation of temperature in 
Casimir effect experiments is almost an impossibility, es¬ 
pecially at low temperatures where vibrational noise is 
difficult to remove due to the boiling of cryogenic liquids 





















5 


(^ . an experimental test of the effects of weak local¬ 
ization on the Casimir effect could more easily be per¬ 
formed at fixed low temperature with varying magnetic 
field, looking for the effect shown in Fig. [5] Only weak 
magnetic fields would be necessary for such an experi¬ 
ment, as applying a magnetic field as weak as tens of 
gauss perpendicular to the plates would be enough to re¬ 
duce the effect by a significant percentage. We propose 
that an experimental test of these effects could be per¬ 
formed in the normal plate-sphere geometry with a very 
thin metallic film at a fixed separation, at a fixed low 
temperature, and with a varying weak magnetic field. 
While the exact numerical values the forces measured in 
this geometry are almost guaranteed to differ from the 
results we find, the general trends in the temperature 
and magnetic field dependence of the force are expected 
to remain. 


V. MESOSCOPIC DISORDER FLUCTUATIONS 

When considering disordered systems one must deter¬ 
mine when to perform averaging over disorder potential 
realizations. Different realizations of the disorder poten¬ 
tial will give different Casimir energies, and local fluc¬ 
tuations in the disorder will cause certain patches on 
each plate to vary in how attractive they are-very similar 
to t he p henomenon of universal conductance fluctuations 
(UCF)-leading to a self-averaging of the Casimir 
energy between two macroscopic plates. This argument 
would imply that instead of carrying out the averaging 
procedure on the linear response 11, which gives the Lif- 
shitz formula with the Drude model, we should perform 
averaging over the entire Casimir energy itself. Indeed, it 
is well known that using the Drude model in the Lifshitz 
formula can not be entirely correct, as it leads to a vio¬ 
lation of the Nernst heat theorem [ij, [lJ|, finding a non¬ 
zero entropy in the limit of zero temperature. In practice, 
however, it is not possible to consider an exact disorder 
potential or to perform the averaging procedure over the 
entire Lifshitz formula, and it is unknown if the simplifi¬ 
cation of using the disorder averaged linear response (i.e. 
the Drude model) in the Lifshitz formula is a legitimate 
approximation for Casimir force measurements despite 
the entropy problem. Another way of phrasing this issue 


£c[Ili,U2]--kBTj2' [t[dn 

{ujn} x=te: 


is that the approximation 

{£c[U])=£4{n)]+6£c^£c[{n)] (9) 

certainly violates the Nernst theorem, but it is unclear if 
it nonetheless closely approximates the exact expression 
for the Casimir energy one would obtain if disorder were 
to be treated exactly or if averaging were done at the 
appropriate stage of the calculation. 

Here we calculate what effect fluctuations from the av¬ 
erage in any particular realization of a disorder potential 
have on the Casimir energy at low temperature, where 
conductance fluctuations are strongest. We start from 
a microscopic version of the Lifshitz formula in position 
space, 

f,[ni,n 2 ] = fcsT^ 'Trln(l-M), (10) 

where 

M =J dridr2dr3lii{r,ri)D{ri,r2)U2{r2,r'^)Dir3,r'). 

Here, D is the photon propagator, which in M con¬ 
nect the screened response of one plate to the other, 

and Hi is the RPA screened electromagnetic linear re¬ 
sponse functions for plate f, schematically given by 
(1 — IIili(0))“^IIi, where H is the unscreened linear re¬ 
sponse function and D{0) is the photon propagator along 
the plate. From this point on we will consider plate 1 to 
be disordered with a particular disorder realization and 
for simplicity we will assume that plate 2 is homogeneous 
and not disordered. We will further take both plates to 
be two-dimensional. For two-dimensional plates, the lin¬ 
ear response function and photon propagator are 2x2 
matrices, with the components of H being proportional 
to the ac conductivity of the plates. The trace, Tr, is a 
generalized trace over both this matrix structure and the 
position labels of the function M. 

We are interested in the case of metallic plates without 
a Hall effect so both response functions are proportional 
to the identity matrix. The photon propagator is diago¬ 
nal as well. The matrix trace in Eq. m becomes trivial, 
leaving us with a sum over photon polarizations. We 
make the further approximation that the Casimir energy 
is well described by just the first term obtained from ex¬ 
panding the logarithm. 


nf (r4, ri)D^{ri,r2)U2ir2, r3)D^{r3, r4), (II) 


where now we label the two photon polarization with tion can be represented as in Fig. [S] 

the superscript X. Diagrammatically, this approxima- r n • j i x- n 

° ensure the following procedures are analytically 
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FIG. 5. The lowest order approximation to the Casimir energy 
given in Eq. m- The grey ovals represent the RPA screened 
linear response functions, and the wavy lines represent photon 
propagators. 


It is related to the size of the fluctuations of the conduc¬ 
tivity, 6a^, (or equivalently in 2D, the conductance) in 



tractable, we must make one further simplifying approx¬ 
imation. We assume that the response function for plate 
1, Hi, which depends on an exact disorder realization, 
can be written as Hi = (Hi) -|- dlli, i.e. the exact re¬ 
sponse function can be written as the disorder averaged 
(Drude) response plus another small term to account for 
the particular disorder realization, here called <5111. (Note 
that this ^IIi is distinct from the function of similar name 
given in Eq. ([S]) and Eq. ([5]) which define a correction to 
the correlation function obtained after disorder averag¬ 
ing.) With this notation, we now expand the RPA on 
plate 1 to first order in i5ni as, 


Hi 


(Hi) 


1- (ni)D 


(ni)D 


6Ui 




= nf 


1 


i-(ni)D(o) 


i-(ni)D (i-(ni)D) 

(i-hnfD(o))(5ni. (12) 


=ri 


This is the form of Hi that is used in Eq. m- 

We now look to the probability distribution of the 
Casimir energy due to fluctuations in the disorder re¬ 
alization in plate 1, now contained entirely within the 
function (5ni, which is given by 


(13) 

,r4)(5ni(r3,r4) 


This expression makes use of the disorder averaged cor¬ 
relator of two unaveraged response functions for the dis¬ 
ordered plate, which can be written as, 

K^^{ri,r2,r3,ri) = ((5ni(r-i,-r2)(5ni(r3, r4)) , (14) 

and is also given diagrammatically in Fig. [B] The func¬ 
tion is very similar to the central object of interest 
considered in the context of universal conductance fluc¬ 
tuations [2ll - l^ . and it is calculated in the same manner. 


Vs^ [£] = {5 {£,-£)) 

= J V{SUi) J 

1 f 

9 / n'^D<5ni(ri,r2)AA( 


X exp 


2=1 


FIG. 6. The primary diagram giving the correlation of dis¬ 
order fluctuations, as defined in Eq. m. The components 
of the diagrams have the same meanings as given in Fig. [S] 
Diagrams containing more diffusons are found either not to 
contribnte to the correlator or are found to have a contribu¬ 
tion 0{lj€FT) smaller. 


related to the conductivity. The only difference between 
the calculation of this function here and in the context 
of UCF is that the latter is primarily concerned with 
conduction of electrons through a system with attached 
leads, usually at zero temperature, while we consider a 
system with no leads at finite temperature. As such, 
most of the qualitative properties of conductance fluctu¬ 
ations apply in our analysis of fluctuations in the Casimir 
energy as well, though the exact form of differs by 
small numerical factors. With this insight, we can al¬ 
ready draw several conclusions about the nature of the 
distribution we will obtain from Eq. (US}. Most impor¬ 
tantly, for weak disorder we can expect fluctuations of 
the Casimir energy around the average value to be small 
since conductance fluctuations are small in good metals: 
i5cr^/cr^ ~ l/(eFT)^- Additionally, we could expect the 
size of the fluctuations to be reduced by a factor of 2 
if a magnetic field were applied to the sample. This is 
because the diagram for ATf ^ given in Fig. [6] gives the 
same contribution at zero magnetic field if all diffusons 
are replaced with cooperons, but the cooperon contribu¬ 
tion is suppressed in magnetic fields in the same way as 
the weak localization correction to the conductivity. 

In order to evaluate Eq. m, we perform a saddle point 
approximation on the functional integral over the disor¬ 
der fluctuation, (^IIi, which after a straightforward cal¬ 
culation gives. 




1 

, exp 

V27rlE2 


{£-£i 


Drude 


21T2 



(15) 


where we have defined the quantities the average, 

and W, the width of the energy distribution. We find 
that the average energy is given by the same expression 
as in Eq. m, but with the substitution IIi —^ Ilf, i.e. 
replacing the exact unaveraged response function IIi with 
the disorder averaged (Drude) response. Therefore, the 
average is simply an approximation of the exact 

Drude result. Additionally, we find that the square of 
the width of the distribution can be written explicitly as. 
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n drrdr [^ (ri, r2)nf (ra, (rg, r4) x 

i^l X,Y 




(16) 



FIG. 7. The diagrams giving the width of the distribution, 
explicitly given in Eq. Ha. 
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which can be represented diagrammatically as in Fig. [3 
The multiple diagrams in this figure result from an ex¬ 
pansion of the r 1 factors of Eq. (HU), 

(i + iiYd^{o))(i + uYd^{o))k-^ 

= (i + n?7^^(o)) (i + n?i:>^(o)) k ^- 

We compute these expressions numerically in the same 
way that we calculate the Casimir pressure in Sec.|TVl For 
both plates we use the Fermi energy and electron mass of 
gold, and we consider plate 1 to be disordered while plate 
2 is a clean plasma plate. We use the material parame¬ 
ters for gold in plate 2 so that in the limit of weakening 
disorder we are left with identical plasma plates. We vary 
the parameter r to determine the dependence of W and 
Eq, and the numerical results for their ratio are fit to 
the expected functional dependence, as shown in Fig. [S] 
In the parameter range we are interested in, we find the 
result. 


W 


cDrude 

^0 


w 


pplasma 



err 


(17) 


In this expression Ci ~ 0.096 is a distance independent 
constant and jg the Casimir energy between two 

clean plasma model plates calculated in the same approx¬ 
imation as given in Eq. (EH)- In the same way, we 

also find this ratio’s dependence on the distance between 
the plates. It suffices to consider only the first term for 
this purpose, since the second term in Eq. dm) has no 
dependence on a. We find. 


W 


pplasma 



(18) 


FIG. 8. (Color online) The fit of numerical data for the quan¬ 
tity to the expected functional dependence given in 

Eq. (I17II . The black dots are the numerical data, the dashed 
blue line is the fit function, and the dotted red line is the 
asymptotic value which has no dependence on r 

in the leading approximation. 


where C 2 ~ 0.038 is another constant independent of 
both r and a. The form of the disorder dependence in 
Eq. (1171) is expected since a weakening of disorder, inter¬ 
preted as an increase of the scattering time r, will make 
a disordered plate more like a plasma plate. Therefore, a 
very large scattering time should give a very good approx¬ 
imation to the plasma result. Note, however, that the 
complete removal of disorder through the limit r —^ 00 
has no physical meaning at this point in the calculation, 
since W has already necessarily been calculated in the 
presence of disorder. We can see from these two ex¬ 
pressions that the distribution will be relatively sharply 
peaked, in the sense that 1, for plates that 

are not too close together and are in the disorder regime 
I/cft "C 1, as we have considered thus far. 

We can get a better understanding of how peaked the 
energy distribution is around its average value by com¬ 
paring its width IF to a smaller relevant energy scale, 
£;Drude _ £-piasma^ combining Eq. (ITTl) and Eq. (fTSl) . 
We obtain. 


IF 


icDrude ^plasma 

\^0 ^0 



(19) 


We see that the nature of the energy distribution Eq. dm 
depends on the two dimensionless quantities hcl{epa) 
and cr/a. Both of these dependencies can be understood 














FIG. 9. (Color online) A plot of the distribution Eq. (1151) given 
for several values of a for a constant value of r = 4.5 x 10“^^ s, 
corresponding to Z = 60 nm. The values of a are 250 (solid 
blue), 400 (dashed green), 800 (dash-dotted yellow), and 
1600nm (dotted red). The average and width W are 

calculated numerically using Eq. (Hill and Eq. (11611 . The plots 
are scaled so the distributions are all the same height, and so 
I^Drude _ £.piasma| ai.^\,a,ys the Same width. Also indicated is 

the value of One sees that as a is increased the distri¬ 

bution becomes sharply peaked even compared to the small 
energy scale set by the difference from the plasma model. 

intuitively. The dependence on hc/^epa) can be under¬ 
stood as arising from the relevant photonic energy scale. 
The most important photons are those with wavelength 
equal to twice the distance between the plates, and when 
this distance is large, these long wavelength photons are 
able to average over larger areas of the plates, reducing 
the effect of local fluctuations. The dependence on cr/a is 
similarly straightforward. It is a comparison of two time 
scales: the impurity scattering time, r, and the time for 
photons to traverse the distance between the plates, a/c. 
When the ratio is small, electrons will have many impu¬ 
rity scattering events before interacting with a photon, 
so any effects due to impurities will be very important. 

There are several regimes we can now explore. Here, 
we will always consider plates of the same material, so 
the Fermi energy is a fixed parameter and we can only 
vary a and r. First, if the plates are very close, meaning 
hc/{eFa) is large, then the distribution is very wide re¬ 
gardless of the size of r. Second, if r is large compared 
to a/c, meaning that photons interact with any given 
electron many times between impurity scattering events, 
then the distribution is again very wide, regardless of the 
size of hc/{eFa). The only regime in which the distribu¬ 
tion is very sharply peaked is when both dimensionless 


parameters are small. This requires that the plates are 
much farther apart than both length scales hc/ep and 
CT, so that each electron undergoes many impurity scat¬ 
tering events between photon interactions and effect of 
disorder is more pronounced, but also so that long wave¬ 
length photons are most important, averaging out the 
disorder fluctuations. 

Ultimately, this result means that for a given level of 
disorder, we can always go to large enough plate separa¬ 
tions so that relatively small local fluctuations in the dis¬ 
order potential of metallic plates are not likely to greatly 
affect the Casimir energy, as shown in Fig. [5] The diffi¬ 
culty here is that for large values of the inelastic scat¬ 
tering time r, the distance at which the distribution 
becomes very sharply peaked may be so large that the 
Casimir effect itself will become unmeasurably small. We 
note now that the values of the parameters r and a used 
in Sec. Hill giving the results in Sec. US give a distribu¬ 
tion very sharply peaked around its average, so we are 
justified in our use of the Drude model despite any of the 
stated concerns over the disorder averaging procedure. 


VI. CONCLUSION 

As we have shown, the weak localization correction 
to the Drude model at low temperatures may give 
the Casimir pressure a nontrivial dependence on both 
temperature and applied perpendicular magnetic field. 
Moreover, we find that, for low enough temperatures, 
WL effects changes the Casimir pressure from the ex¬ 
pected value without WL by an amount greater than the 
difference between the Drude and plasma model predic¬ 
tions. Since these effects are not applicable in a model 
of a 2D plate without disorder, i.e. the plasma model, a 
high precision experimental test measuring this tempera¬ 
ture or magnetic field dependence would give a definitive 
indication of whether a diffusive model truly describes 
the behavior of electrons in Casimir experiments. 

Additionally, we explore the effect that fluctuations in 
the disorder potential can have on the Casimir energy 
and the validity of using the Drude model considering 
that the correct averaging procedure would give a result 
that differs from the Drude model by the inclusion of 
nonlocal disorder fluctuation contributions. We find that 
for a given level of disorder, one can always overcome the 
effects of fluctuations by holding the plates far enough 
apart, which justifies the use of the Drude model. 
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